Local control of entanglement in a spin chain 
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In a ferromagnetic spin chain, the control of the local effective magnetic field allows to manipulate 
the static and dynamical properties of entanglement. In particular, the propagation of quantum 
correlations can be driven to a great extent so as to achieve an entanglement transfer on demand 
toward a selected site. 

PACS numbers: 03.67.Hk, 03.67.Mn, 05.50.+q 

When the superposition principle is applied to an (at 
least) bi-partite system, highly non local and purely 
quantum correlations (entanglement) can appear among 
the parties. This constitutes a crucial resource for 
many applications in quantum communication [1]. Con- 
sequently, the problem of entanglement distribution has 
become of central interest: quantum correlations are gen- 
erated by local interactions; therefore methods are re- 
quired to transfer either the entangled particles or their 
state at a distance. It has been shown theoretically that 
spin chains are efficient quantum channels for short dis- 
tance entanglement transfer |2|; thus the ability to ma- 
nipulate the propagation of entanglement in a spin chain 
can be very important and it has already been shown 
that breaking the translational invariance of the chain 
can produce very interesting results in this respect [1, . 

In this paper, we show that a control of both static and 
dynamical properties of entanglement can be achieved by 
acting locally to modify the level spacing of some qubits. 
In the magnetic language, we analyze a system subject 
to a spatially inhomogeneous magnetic field; that is, a 
system with 'diagonal defects' ('impurities'), 

In contrast to the usual case in which local actions can- 
not affect non-local physical quantities; here, due to the 
spin-spin interaction, the local control of the magnetic 
field modifies both the distribution of correlations in the 
ground state and the entanglement propagation along the 
chain. Indeed, spatial inhomogeneities of the external 
field lead to an Anderson-like localization of entangle- 
ment 0, and this can occur even for a single defect [7], 
giving rise to a mirror-like effect in the entanglement 
propagation [8]. Contrary to the usual description of 
localization phenomena, we do not conceive the impurity 
just as a bit of disorder in the system, but rather intend 
the modification of the local level spacing as a knob to i) 
control the content of "static" (ground state) entangle- 
ment, and a) drive its propagation along the chain. 

We consider a 1-D XX spin-^ closed chain, placed in an 
external magnetic field which is homogeneous everywhere 
but for two defect sites h and This model is described 
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where uoq is the level spacing of each qubit, except for 
those residing at sites Z^, which have level separation 
uoq + a^. Henceforth, we shall take J = 1 and use the 
ferromagnetic coupling constant as our energy unit (we 
have also set h= 1). Furthermore, we set to zero the en- 
ergy of the completely separable eigenstate |0)*^^, which 
is the unperturbed ground state for cjq > 1, and which is 
still an energy eigenstate in the presence of defects. 

This model can be solved exactly via the Jordan- 
Wigner (JW) transformation, which maps the spin chain 
into a spin- less fermion system [9]. However, we do not 
need the general solution here, since we will deal with 
states having at most one tilted spin (in the JW lan- 
guage, states lying in the single particle subspace), [l^. 
In the continuum limit ^ oo, we solve the model (re- 
stricted to the single particle sector of its Hilbert space) 
using the Green operator technique f\. 

To this end, we consider the Green operator Gq de- 
scribing the homogeneous chain. It is known 0, @1 that 
Go displays a branch cut on the real axis in the complex 
energy plane for energies E ^ [uoq — I^uoq ^ 1]. This cut 
signals a continuous energy band, which survives even in 
the presence of the defects. 

Given Go, the full (single particle) Green operator is 



G — Go + GoTGo, 
where the T-matrix is given by 
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where |n) = |0) ~ is the state with one spin 
down (or, equivalently, one fermionic excitation) located 
at site n, and where the scattering coefficients are 
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The T-matrix describes multiple scattering events of the 
single particle excitation at the two defects, and it is 
precisely the re-summation of the various scattering am- 
plitudes which gives rise to the denominator in Eq. (|4]). 
The existence of zeros of this denominator is crucial since 
it implies that, besides the energy band discussed above, 
the model with defects displays some (at most two, in 
fact) discrete energy levels. 

They lay above or below the energy band depending 
on the values of the defect fields, and their eigen-energies 
can be obtained analytically in some special cases. For 
example, for nearest neighbors defects, with distance d = 
1^2 - I = 1, by setting xioc = Eioc - one gets 
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Only one of these two solutions (the one with the plus 
sign in front of the square root) exists if the defect 
strengths satisfy the relation + > 1. That is, 
the parameter space is divided in two regions, character- 
ized by the number of discrete energy eigenvalues (which 
can be one or two). 

If we restrict ourselves to the ordered phase of the un- 
perturbed chain, uoq> 1^ and if we consider only ai > 0, 
the lowest (or the only existing) localized level becomes 
the fundamental state of the perturbed problem. 

Many of these features obtained for nearest neighbors 
defects are generic and do not depend on their distance. 
From the analytic properties of the Green operator, 0, 
one can show that i) at most two localized states are 
present, whose position with respect to the energy band 
depends on the sign of the a's, ii) a region in the ai- 
a2 plane exists in which only one discrete state is found. 
However, this one-eigenstate region becomes thinner and 
thinner as d increases, iii) the localized states always 
have the structure of single excitation states of the form 
li^ioc) = X^n^^l^)- amplitudes 6n, obtained from 

the residues of G, can be expressed in terms of an in- 
verse localization length ^ = — \n.[—xioc — V^foc 
a function "bi-localized" around the two defects: 

6n = const(Kie-^l^-^^l +K2e-^l^-^^l). 
The coefficients Ki are given by 
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Their ratio indicates the relative weights of the two local- 
ization region. For ai ^ a2 one gets \K2 \ <C \Ki \ and the 
discrete levels are localized around /i, while for ai <C a2 
the localization center is I2 since \K2\ ^ For equal 

defect strengths, ai =0^2, one finds K2 = where 
the upper (lower) sign refers to the lower (higher) of the 
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FIG. 1: (Color online) Ground state concurrence between the 
/i-th and the n-th spin of the chain. The entanglement can 
be remotely controlled by changing the local field at the other 
defect. The line with open circles (boxes) corresponds to 
tti = tt2 = 1-5 {ai = 2,0^2 = 1-5). For equal defect strengths, 
spin h is entangled both with its own neighborhood and with 
the other defect's one. In the asymmetric case, quantum cor- 
relations only survive within the localization region. The inset 
shows the concurrence between the two defects as a function 
of their distance for various values of defect fields. From below 
to above, the plots correspond to ai = a2 = 0.25; 0.5; 1; 2. 



two eigenstates. Thus, if ai =0^2, the two discrete states 
are given by equal- weight coherent superpositions of two 
localized parts centered on the two impurities. These 
states are highly entangled and display strong quantum 
correlations between the defects and their neighborhoods. 

Even for quite small values of the defect fields, the 
localization length is smaller than the inter-site spacing 
(already ai = 0^2 = 1 gives < 1, for any distance 
d = \li — l2\). This implies that the two discrete eigen- 
states are approximately given by the Bell combinations 
(|01) ± 1 10)), with the rest of the chain almost completely 
factorized in the state |0). It is noteworthy that this 
structure does not depend on the distance between the 
two defects, see the inset of Fig. (TJ where the concur- 
rence Ci^i^ = 2 1 61 62 1 between the defects is shown as a 
function of d = |/i — /2|- Thus, a long distance bi-partite 
entanglement can be obtained in the ground state of the 
chain. A similar behavior has been found in Ref. (Tlj , 
with the difference that in our case this is a bulk property 
rather than a surface effect. 

For generic values of the defect amplitudes, these states 
display entanglement between any pair of spins residing 
near each of the two defects, with the peculiarity that 
the pairwise entanglement for two spins around the same 
defect depends on the value of the local magnetic field at 
the other defect. This is illustrated in figure ([1]), where it 
is shown that the entanglement around li is modified by 
changing the strength of local field at ^2, thus achieving 



3 




C,. 



r 



(XlO')t 




(X 10') t 



FIG. 2: (color online) Rabi oscillations of the entanglement 
between the two defects for the case ai = 0^2 = 1-5. 



a remote entanglement control. 

The bi-local character of the discrete levels strongly 
affects the transport of entanglement along the chain. In 
particular, we consider the possibility of using the chain 
to send one partner of a maximally entangled pair. We 
assume that the spin at the sender site s is prepared in 
a singlet state with an external (un-coupled) qubit. The 
interaction between the spins causes a transfer of entan- 
glement along the chain. Ideally, after a given transmis- 
sion time t, one would like to get a singlet between the 
external qubit and the one residing at a receiving site r. 
To characterize the quality of the transmission, we evalu- 
ate the (final) concurrence between the external and the 
r-th qubits, denoted by Cr{t). 

This is given by Cr{t) = \fs^r{t)\, [2], where fs^r(t) 
is the amplitude for the transfer of a fermionic excitation 
from site s to r. This is expressed in terms of the retarded 
Green operator as 



fs^r = ^ < 
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The first term describes transport mediated by localized 
states, while the second one gives a spin- wave mediated 
transfer, [12]. We analyze them separately. 

The first contribution is effective only within a region 
of length around the two defect sites. This gives the 
noteworthy possibility of transmitting from one neigh- 
borhood to the other. This effects is illustrated in Fig. 
(j2j), where the entanglement is shown to jump from one 
defect to the other. Indeed, if the sender site coincides 
with one of the defects (say /i), than the system evolves 
by performing almost perfect Rabi oscillations between 
the two discrete localized states, with a Rabi frequency 
given by uor = {Eiod — E10C2)' This can be understood 
by noticing that the initial singlet state between li and 
the external qubit can be approximately written as 
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FIG. 3: (color online) Entanglement bouncing between the 
two defects which act like (non-perfect) mirrors. The location 
of the impurities is indicated by the two arrows and the defect- 
strengths are ai = a2 = 1.5. 



which implies that the concurrence between the li-th 
qubit and the external one is an harmonic function: 



Ci^{t) \ cos ujRt\, Ci^{t) \ sin ujRt\. 



(7) 



The second term in Eq. (|6]) is an integral over the energy 
band, parameterized as E = coq — cos 0. It contains the 
state amplitudes of the continuous energy band, which 
can be written in terms of the retarded Green and T 
operators: 

gn{E) = {n\ [l+G+T+j \ME)) . 

with the un-perturbed states such that (n|?/^o(^)) = 
e*^^/>/27r. These states represent distorted spin waves of 
the system. They are the stationary scattering states of 
the single-particle Hamiltonian and can be constructed, 
starting from the usual magnon excitations, by includ- 
ing corrections due to multiple scattering at the defects. 
For moderate values of the a's, the distortion is not 
that big and the unperturbed plane wave nature can still 
be recognized. As a result, the energy eigenstates form 
(approximate) standing waves between the two defects. 
Their pattern is reflected in the entanglement propaga- 
tion shown in Fig. (|3]) where the sender site is located 
between the two impurities which act as potential bar- 
riers for the spin waves, so that entanglement bounces 
back and forth between these two mirrors. The extreme 
situations is reached when the defects are next to nearest 
neighbors to each other, thus realizing an entanglement 
trap, see Fig. (|4]). Since the mirrors are not perfect, 
[8[, the trapped entanglement decreases with time, the 
superimposed time oscillations being due to the e~*^'°^^ 
factor of the symmetric discrete eigenstate (the only one 
that matters, in this case). Once these oscillations are 
subtracted, the short time behavior of the concurrence is 
found to be parabolic, with a convexity that decreases as 
the defect amplitudes are increased. A long-time (resid- 
ual) trapped entanglement is also present, which is due 
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FIG. 4: (color online) Left: Entanglement trapping between 
the two defects in the case ai = a2 = 2. Right: Time 
dependence of the trapped concurrence for ai — cx2 — 
0.5, 1, 1.5, 2.5, 5 (from below to above). 



to the tails of the localized state. It diminishes with in- 
creasing a's as the localization length does, and becomes 
negligible if is much smaller than the site spacing. 

The Rabi oscillations illustrated above give a mean to 
transfer the entanglement reversibly between the defects. 
This method works independently of their distance, but 
has the drawback that the Rabi period increases with 
distance as T ~ (this can be derived by perturbation 
theory for large a's, see Ref. [13]). 

The form of the discrete states, however, suggests an- 
other, more effective method to achieve entanglement 
transfer on demand between the two defects; namely, the 
adiabatic passage (see [14] for a related proposal in which 
the coupling strength is varied instead of the local field). 

The idea is to change the defect strengths adiabatically, 
so that the ground state of the system, having the general 
form given in Eq. ([5j), is adiabatically changed from the 
initial state \i\)gs(i)) — localized at /i, to the final 
state |'0^s(/)) — 1/2), localized at 12- This can be done 
by modifying the defect fields from the initial values 



If this is done adiabatically, the system always remains in 
its instantaneous ground state, thus realizing the desired 
entanglement transfer provided the initial singlet state 
involves the external and the /i-th qubit. To ensure the 
adiabaticity, the rates of change of the a's have to bo 
much smaller than the difference between the energies of 
the two discrete levels, £^21 = ^ioc2 — ^loci- The most 
dangerous point in this respect (i.e., the smallest £^21) 
occurs at the crossing, when Oi\(i) = 0^2 (^). But, since 
at this point £^21 ^ <^^, if the adiabatic procedure is 
designed such that the crossing occurs for a very small a, 
than the adiabaticity can be preserved even for transfer 
times much smaller than the Rabi period. This procedure 
has the additional advantage of effectively decoupling the 



receiving site from the rest of the chain after the transfer 
has been performed due to its large final local field. 

To summarize, we have discussed how to manipulate 
a spin chain with local control fields, showing that the 
static entanglement can be remotely controlled, and that 
entanglement propagation can be adjusted to a large ex- 
tent in order to achieve transfer on demand. One possi- 
bility to implement this model in a realistic set up is to 
use the method proposed in Ref. [lH to "engineer" spin 
chains with atoms in an optical lattice. The addition of 
external static local electric or magnetic fields should en- 
able the control of the qubit energy level spacing which 
is essential to test our proposal. 

This work has been supported by MIUR through the 
PRIN2005029421 project. 
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